Usenet 1993 July/InfoMagic USENET CD-ROM July 1993.ISO/sources/misc/volume9/draw_jt/part02

Text File | 1989-12-12 | 51.3 KB | 2,135 lines

Newsgroups: comp.sources.misc
From: [e-mail] (Brandon S. Allbery - comp.sources.misc)
Subject: v09i062: draw, part 02/02
Reply-To: [e-mail] (Juha Takala)

Posting-number: Volume 9, Issue 62
Submitted-by: [e-mail] (Juha Takala)
Archive-name: draw_jt/part02

[This is "dr_lib", a library to produce plots in various formats including
plot(3) and HPGL formats, and "draw", a simple front-end for the library.
See the README for more information.  ++bsa]

#! /bin/sh
# This is a shell archive.  Remove anything before this line, then unpack
# it by saving it into a file and typing "sh file".  To overwrite existing
# files, type "sh file -c".  You can also feed this as standard input via
# unshar, or by typing "sh <file", e.g..  If this archive is complete, you
# will see the following message at the end:
#        "End of archive 2 (of 2)."
# Contents:  dr_lib.c draw.c
# Wrapped by allbery@uunet on Tue Dec 12 21:15:58 1989
PATH=/bin:/usr/bin:/usr/ucb ; export PATH
if test -f 'dr_lib.c' -a "${1}" != "-c" ; then
  echo shar: Will not clobber existing file \"'dr_lib.c'\"
else
echo shar: Extracting \"'dr_lib.c'\" \(31312 characters\)
sed "s/^X//" >'dr_lib.c' <<'END_OF_FILE'
X/*
X
XCopyright (C) 1988, 1989 by Juha Takala, [e-mail]
X
X     This program is free software; you can redistribute it and/or modify
X     it under the terms of the GNU General Public License as published by
X     the Free Software Foundation; version 1.
X
X     This program is distributed in the hope that it will be useful,
X     but WITHOUT ANY WARRANTY; without even the implied warranty of
X     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
X     the file "License" for more details
X
X     You should have received a copy of the GNU General Public License
X     along with this program; if not, write to the Free Software
X     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
X
X
X */
X
X
X/* Description:
X *  dr_lib.c provides interface between device level and my draw(1).
X *  Supported devices are:
X *  - hpgl output for hp7475 plotter
X *  - gl_plot(3) calls for EGA or matrix printer output
X *  - plot(3) interface
X *
X */
Xstatic char *rcsid =
X  "$Id: dr_lib.c,v 2.4 89/12/08 14:00:34 jta Exp Locker: jta $";
X
X/* define this, if your plotter connection is 'y-cable' (not direct) */
X/* #define YCABLE            /*  */
X
X/* definede this if you want to ask paper size from plotter */
X/* (should not be used, if you are going to use this package as */
X/* part of pipelines) */
X/* #define ASK            /*  */
X
X#include <stdio.h>
X#include <math.h>
X#ifdef MSDOS
X#include <bios.h>
X#endif
X#include "dr_lib.h"
X
X#define  ESC 0x1b        /* ascii <ESC> */
X
Xstatic int xscale(), yscale();
Xstatic void sends(), invdev();
X
X#ifdef ASK
Xstatic void takes();
X#endif /* ASK */
X
X/*
X * internal global variable declarations
X *
X */
X                /* coordinate transformations */
Xstatic float xuserincr,   yuserincr;
Xstatic float xuserlow,    yuserlow;
Xstatic float xuserhigh,   yuserhigh;
Xstatic float xdevicelow,  ydevicelow;
Xstatic float xdevicehigh, ydevicehigh;
Xstatic float xrelarea,    yrelarea;
Xstatic float xfactor,     yfactor;
Xstatic float shape_factor, size_factor;    /* text 'fatness' and size factors */
Xstatic int lastx=0, lasty=0;    /* device coordinates for last point */
X
X
Xstatic int device = INVDEV;    /* what output device we are working for */
Xstatic int old_linetype = -1;    /* remember last used linetype */
Xstatic int old_pen = -1;    /* ...and pen */
X
Xstatic int errcount, warncount, errfile;
Xstatic char *errfilename = "draw.err";
X
Xstatic FILE *errfilefp;
X
X
X
X/* ********************************************************** */
X/*                                                            */
X/* Following functions are ment to be subroutines to be       */
X/* called from a user supplied main program                   */
X/*                                                            */
X/* ********************************************************** */
X
X/*
X * dr_set_cset()   - select cset
X *
X */
Xvoid dr_set_cset(n)
Xint n;
X{
X    char s[50];
X
X    switch (device)
X      {
X      case HP7475:
X      if (n < 0 || n == 5 || (n > 9 && n < 30) || n > 39) {
X          (void)sprintf(s, "Invalid cset code %d, using 0", n);
X          dr_errmsg(s);
X          n = 0;
X      }
X      (void)sprintf(s, "CS%d;\n", n);
X      sends(s);
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X/*      (void)sprintf(s, "cset not yet for SCREEN");
X      dr_errmsg(s);        /*  */
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X/*      (void)sprintf(s, "cset not valid for PLOT");
X      dr_errmsg(s);        /*  */
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X}
X
X/*
X * dr_set_pen()   - select pen
X *
X */
Xvoid dr_set_pen(n)
Xint n;                /* pen number, 1..6, 0=='put it away' */
X{
X    char s[50];
X
X    switch (device)
X      {
X      case HP7475:
X      if (n < 0 || n > 6) {
X          (void)sprintf(s, "Invalid pen code %d, using 1", n);
X          dr_errmsg(s);
X          n = 1;
X      }
X      if (n != old_pen) {    /* optimize for hpgl... */
X          old_pen = n;
X          (void)sprintf(s, "SP%d;", n);
X          sends(s);
X      }
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X/*      (void)sprintf(s, "dr_set_pen(colour) not yet for SCREEN\n");
X      dr_errmsg(s);        /*  */
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X/*      (void)sprintf(s, "dr_set_pen(colour) not valid for PLOT\n");
X      dr_errmsg(s);        /*  */
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X}
X
X/*
X * dr_goto()  - move pen to point x, y
X *
X */
Xvoid dr_goto(x, y, penup)
Xfloat x, y;            /* where he wants to go */
Xint penup;            /* nonzero == lift pen before moving */
X{
X    char s[60];
X    lastx = xscale(x);
X    lasty = yscale(y);
X    switch (device)
X      {
X      case HP7475:
X      if (penup) sends ("PU");
X      else sends ("PA;PD");
X
X      /* go to point, leave pen as told
X       */
X      (void)sprintf(s, "%d,%d;\n", lastx, lasty);
X      sends(s);
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      if (penup)    n_movepen(lastx, lasty);
X      else        n_draw(lastx, lasty);
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      if (penup)    move(lastx, lasty);
X      else        cont(lastx, lasty);
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X    return;
X}
X
X
X/*
X * dr_draw_line()  - draw line (x1,y1) .. (x2,y2)
X *
X * connect the two data points with line
X */
Xvoid dr_draw_line(x1, y1, x2, y2)
Xfloat x1, y1, x2, y2;        /* data points */
X{
X    char s[60];
X
X    lastx = xscale(x2);
X    lasty = yscale(y2);
X    switch (device)
X      {
X      case HP7475:
X      /* lift pen,
X       * go to first point,
X       * drop pen there,
X       */
X      (void)sprintf(s, "PA;PU%d,%d;PD;", xscale(x1), yscale(y1));
X      sends(s);
X
X      /* draw line segment, lift pen
X       */
X      (void)sprintf(s, "PD%d,%d;PU;\n", lastx, lasty);
X      sends(s);
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      n_line(xscale(x1), yscale(y1), lastx, lasty);
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      line(xscale(x1), yscale(y1), lastx, lasty);
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X    return;
X}
X
X
X/*
X * dr_draw_array()  - draw y = f(x)
X *
X * connect data points with line
X */
Xvoid dr_draw_array(x, y, n)
Xfloat x[], y[];            /* data points */
Xint n;                /* number of points in x[] and y[] */
X{
X    int i;
X    char s[60];
X
X    /* if we really have something to draw
X     */
X    if (n < 1)
X      return;
X
X    switch (device)
X      {
X      case HP7475:
X      /* lift pen,
X       * go to first point,
X       * drop pen there,
X       */
X      (void)sprintf(s, "PA;PU%d,%d;PD;", xscale(x[0]), yscale(y[0]));
X      sends(s);
X
X      /* draw line segments
X       */
X      for (i=1; i<n; i++) {
X          if (i % 4 == 0) sends("\n");
X          (void)sprintf(s, "PD%d,%d;", xscale(x[i]), yscale(y[i]));
X          sends(s);
X      }
X
X      /* finally lift pen
X       */
X      sends("PU;\n");
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      n_movepen(xscale(x[0]), yscale(y[0]));
X      for (i=1; i<n; i++)
X        n_draw(xscale(x[i]), yscale(y[i]));
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      move(xscale(x[0]), yscale(y[0]));
X      for (i=1; i<n; i++)
X        cont(xscale(x[i]), yscale(y[i]));
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X    lastx = xscale(x[i-1]);
X    lasty = yscale(y[i-1]);
X    return;
X}
X
X
X/*
X * dr_put_mark()  - make a mark at (x,y)
X *
X * put mark at data point
X */
Xvoid dr_put_mark(x, y, mark)
Xfloat x, y;            /* data point */
Xchar mark;
X{
X    char s[60];
X
X    lastx = xscale(x);
X    lasty = yscale(y);
X    switch (device)
X      {
X      case HP7475:
X      /* if we really have something to plot
X       */
X
X      /* lift pen,
X       * set symbol mode,
X       */
X      (void)sprintf(s, "PA;PU;SM%c;", mark);
X      sends(s);
X
X      /* put symbol at all data points
X       */
X      (void)sprintf(s, "PU%d,%d;", lastx, lasty);
X      sends(s);
X
X      /* reset symbol mode
X       */
X      sends("SM;\n");
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      n_movepen(lastx, lasty);
X      n_grafchar(mark);
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      move(lastx, lasty);
X      *s = mark;
X      *(s+1) = '\0';
X      label(s);
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X    return;
X}
X
X
Xvoid dr_set_shape(x)
Xfloat x;
X{
X    if (x > 0.0)
X      shape_factor = x * xrelarea/yrelarea;
X    return;
X}
X
X
Xvoid dr_set_size(x)
Xfloat x;
X{
X    if (x > 0.0)
X      size_factor = x;
X    return;
X}
X
X
X/*
X * dr_put_text()  - put some text
X * plot(3) only has one size and orientation, what can you do!
X *
X */
Xvoid dr_put_text(x, y, xofs, yofs, angle, size, s)
Xfloat x, y;            /* where to put it */
Xfloat xofs, yofs;        /* offset of start position, measured in */
X                /* character size units */
Xfloat angle;            /* radians, 0.0 == horizontal */
Xfloat size;            /* relative size of text */
Xchar *s;            /* the text */
X{
X    char ss[60];
X    int senddir;
X    char dir[30];
X
X    if (size <= 0.0) {
X    sprintf(ss, "invalid size: %f");
X    dr_errmsg(ss);
X    return;
X    }
X    size *= size_factor;
X
X    switch (device)
X      {
X      case HP7475:
X      senddir = (angle != 0.0);
X      if (senddir) {
X          if (angle == M_PI_2)
X        sprintf(dir,"DI0,1");
X          else if (angle == -M_PI_2)
X        sprintf(dir,"DI0,-1");
X          else
X        sprintf(dir,"DI1,%f", tan(angle));
X      }
X
X      /* lift pen, go to specified position, set direction & size,
X       */
X
X      (void)sprintf(ss, "PA;PU%d,%d;%s;SR%4.2f,%4.2f;\n",
X            xscale(x), yscale(y),
X            (senddir ? dir : ""),
X            size * 0.75 * shape_factor * xrelarea,
X            size * 1.5 * yrelarea);
X      sends(ss);
X
X      /* go to offset,
X       * put text,
X       * restore 0 degree direction.
X       */
X      if ((xofs != 0.0) || (yofs != 0.0)) {
X          (void)sprintf(ss, "CP%4.2f,%4.2f;\n", xofs, yofs);
X          sends(ss);
X      }
X      sends("LB"); sends(s); sends("\03\n");
X
X      /* correction for y-offset, if necessary, to leave pen on same */
X      /* height but after the text */
X      if (yofs != 0.0) {
X          (void)sprintf(ss, "CP0,%4.2f;\n", -yofs);
X          sends(ss);
X      }
X
X      if (senddir)
X        sends("DI1,0;\n");
X/*      lastx = ....;
X      lasty = ....;        /* hp7475 does not really need these, */
X                /* it remembers the pen position */
X
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      n_movepen(xscale(x), yscale(y));
X
X#define GL_MODE_255_HAS_BUG
X#ifdef GL_MODE_255_HAS_BUG
X      /* DT Lewis!!  Wake up! */
X      /* For some reason gl_lib works strangely with GL_MODE == 255 */
X      /* (the output is  (almost) microscopic small and weird shape). */
X      /* We can compensate that strangeness here by modifying size */
X      /* and sape factors for g_fontctl()...  The values 3.0 and 0.4 */
X      /* are experimental results (by [e-mail]). */
X      if (strncmp(getenv("GLMODE"),"256",3)==0)
X        g_fontctl(3.0*size*yrelarea, 0.4*shape_factor, 1.0, angle, 0.0);
X      else
X#endif
X#undef GL_MODE_255_HAS_BUG
X        g_fontctl(size*yrelarea, shape_factor, 1.0, angle, 0.0);
X      yofs = -0.5 - yofs;
X      n_ch_ofs(xofs, yofs);
X      n_grafstr(s);
X      n_ch_ofs(0.0, -yofs);
X      g_fontctl(size, 1.0, 1.0, 0.0, 0.0);
X/*      lastx = ....;
X      lasty = ....;        /* this does not really need these, */
X                /* it remembers the pen position */
X
X      break;
X#endif
X#ifdef PLOT
X#define PLOT_CH_XSIZE 70
X#define PLOT_CH_YSIZE 100
X      case PLOTLIB:
X      {
X          int i;
X          float xsize, ysize;
X          float xbase, ybase;
X          float xincr, yincr;
X          float xofs1, yofs1;
X
X          xsize = PLOT_CH_XSIZE * size * xrelarea;
X          ysize = PLOT_CH_YSIZE * size * yrelarea;
X
X          yofs += 0.4;    /* middle --> botom */
X          xofs1 = (xofs * xsize * cos(angle) * shape_factor -
X               yofs * ysize * sin(angle));
X          yofs1 = (xofs * xsize * sin(angle) * shape_factor +
X               yofs * ysize * cos(angle));
X          xbase = xscale(x) + xofs1;
X          ybase = yscale(y) + yofs1;
X
X          xincr = cos(angle) * xsize * shape_factor;
X          yincr = sin(angle) * ysize;
X
X          /* This is how we simulate different printing directions */
X          /* and character sizes: by moving to the position of every */
X          /* single character one at a time */
X          *(ss+1) = '\0';
X          for (i=0; *ss=*s; s++,i) {
X          move((int) (0.5 + xbase + i * xincr),
X               (int) (0.5 + ybase + i * yincr));
X          label(ss);
X          i++;
X          }
X
X          lastx = (int) (0.5 + xbase + i * xincr - xofs1);
X          lasty = (int) (0.5 + ybase + i * yincr - yofs1);
X          move(lastx, lasty);
X      }
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X    return;
X}
X
X/*
X * dr_small_circle()  - draw small circle at current point.
X * This routine needs the help of lastx and lasty in case of plot(3) library.
X *
X */
Xvoid dr_draw_circle()
X{
X    switch (device)
X      {
X      case HP7475:
X      sends("CI35;\n");      break;
X#ifdef GL_LIB
X      case SCREEN:
X      n_ellipse1(100, 100);      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      circle(lastx, lasty, 20);      break;
X#endif
X      default:
X      invdev(); break;
X      }
X}
X
X
X/*
X * dr_finish()  - terminate the plotting
X *
X * This routine must be called after all plotting activity is done
X * to restore the Y-cable connection so that terminal gets control
X */
Xint dr_finish(clean_scr)
Xint clean_scr;            /* != 0 if we have clear screen */
X{
X    char s[20];
X
X    switch (device)
X      {
X      case HP7475:
X      /* deselct pen,
X       * reset Y-cable connection
X       *   (<esc>.Z is ignored if not Y-cable connecton )
X       */
X      (void)sprintf(s, "SP;%c.Z\n", ESC);
X      sends(s);
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      if (clean_scr && (strncmp(getenv("GLMODE"), "256", 3) != 0)) {
X#ifdef UNIX
X          FILE *fp = fopen("/dev/tty", "r");
X          fgets(s, 2, fp);    /* wait for user response */
X          fclose(fp);
X#endif /* UNIX */
X#ifdef MSDOS
X          _bios_keybrd(_KEYBRD_READ); /* wait for keybrd. input forever */
X#endif /* MSDOS */
X      }
X      {
X          int status;
X          status = g_finish();    /* do clean screen */
X          if (status != 0)
X        dr_errmsg("dr_lib(): problems with g_finish().");
X      }
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      move(0,0);
X      closepl();
X      break;
X#endif
X      default:
X      break;        /* call to invdev() --> recursion */
X      }
X
X    /* check errors, tell if any or remove error log file
X     */
X    if (errcount == 0 && warncount == 0)
X      (void)unlink(errfilename);
X    else {
X    if (errcount)
X      (void)fprintf(stderr, "%d error%s,",
X            errcount, (errcount == 1 ? "" : "s"));
X    if (warncount)
X      (void)fprintf(stderr, "%d warning%s,",
X            warncount, (warncount == 1 ? "" : "s"));
X    (void)fprintf (stderr, " see file \"%s\".\n", errfilename);
X    }
X    /* return something useful
X     */
X    return (errcount);
X}
X
X
X/*
X * dr_start  - init the plotting device
X *
X * This routine must be called before any plotting activity is done
X * to set up the Y-cable connction so that plotter starts processing
X * or in case of direct connection, open the device
X */
Xvoid dr_start(xlow, ylow, xhigh, yhigh, dev)
Xint dev;            /* plotting device & paper size */
Xfloat xlow, xhigh;        /* drawing area, relative to.. */
Xfloat ylow, yhigh;        /* ..whole paper, that means: */
X                /* 0.0, 0.5, 0.5, 1.0 will use */
X                /* upper left quarter of paper */
X{
X    int p1x, p1y, p2x, p2y;
X    float t1, t2;
X    char s[200];
X
X    /* open error log file
X     */
X    errcount = 0;
X    warncount = 0;
X    if ((errfilefp = fopen(errfilename, "w")) == NULL) {
X    (void)fprintf(stderr, "\nCan't open %s.  Aborting.\n", errfilename);
X    exit(1);
X    }
X    errfile = fileno(errfilefp);
X    xrelarea = xhigh - xlow;
X    yrelarea = yhigh - ylow;
X    shape_factor = yrelarea/xrelarea; /* back to normal 'fatness' */
X    size_factor = 1.0;        /* actual size normalized to area */
X
X    switch (device = (dev & DEVMASK))
X      {
X      case HP7475:
X      /* set Y-cable connection,
X       *   (<esc>.Y is ignored if not Y-cable connecton )
X       * init the plotter,
X       * set handshake mode x-on/x-off,
X       * select pen #1
X       */
X      (void)sprintf(s, "\n%c.Y;IN;", ESC); /* Y-connection, plotter on */
X      sends(s);
X      (void)sprintf(s, "%c.I80;;17:", ESC);    /* enable x-on/x-off x-on=^q */
X      sends(s);
X      (void)sprintf(s, "%c.N;19:\n", ESC);    /* x-off=^s */
X      sends(s);
X#ifdef ASK
X      /* ask picture area from plotter
X       * set up limits (leave some room for axis labels etc.) */
X      sends("OP;");
X      takes(s);
X#else /* ASK */
X       /* we would get one of this */
X      switch (dev & SIZEMASK)
X        {
X        case SIZE_A:
X        (void)sprintf (s, "250,596,10250,7796\n");
X        sends("PS4;\n");
X        break;
X        case SIZE_A4:
X        (void)sprintf (s, "603,521,10603,7721\n");
X        sends("PS4;\n");
X        break;
X        case SIZE_B:
X        (void)sprintf (s, "522,259,15722,10259\n");
X        sends("PS3;\n");
X        break;
X        case SIZE_A3:
X        (void)sprintf (s, "170,602,15370,10602\n");
X        sends("PS3;\n");
X        break;
X        default:
X        dr_errmsg("Invalid paper size.  Aborting.");
X        (void)dr_finish();
X        exit(1);
X        }
X#endif /* ASK */
X      if ( sscanf(s, "%d,%d,%d,%d", &p1x, &p1y, &p2x, &p2y) == 4 ) {
X                /* some room for labels */
X          t1 = p1x + 0.05*(p2x-p1x); /* y-label */
X          t2 = p2x - 0.001*(p2x-p1x); /* right margin */
X          xdevicelow  = t1 + xlow * (t2 - t1);
X          xdevicehigh = t1 + xhigh * (t2 - t1);
X          t1 = p1y + 0.05*(p2y-p1y); /* x-label */
X          t2 = p2y - 0.07*(p2y-p1y); /* name */
X          ydevicelow  = t1 + ylow * (t2 - t1);
X          ydevicehigh = t1 + yhigh * (t2 - t1);
X      }
X      else {
X          dr_errmsg("Can't receive plotter dimensions.  Aborting.");
X          (void) dr_finish();
X          exit(1);
X      }
X
X      dr_set_pen(1);        /* select some pen */
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      p1x = 0;    p1y = 32767;
X      p2x = 32767;    p2y = 0;
X      t1 = p1x + 0.06*(p2x-p1x); /* room y-axis label */
X      t2 = p2x - 0.001*(p2x-p1x); /* no room at right margin */
X      xdevicelow  = t1 + xlow * (t2 - t1);
X      xdevicehigh = t1 + xhigh * (t2 - t1);
X      t1 = p1y + 0.10*(p2y-p1y); /* room for x-axis label */
X      t2 = p2y - 0.05*(p2y-p1y); /* room fo picture's name */
X      ydevicelow  = t1 + ylow * (t2 - t1);
X      ydevicehigh = t1 + yhigh * (t2 - t1);
X
X#define ENV_MODE 0        /* from gl.h; it means that g_init() will */
X                /* look operating mode in environment */
X                /* variable GLMODE */
X      g_init(ENV_MODE);
X#undef ENV_MODE
X      /* dr_set_pen(1);        /* select some pen */
X      break;
X#endif
X
X#ifdef PLOT
X      case PLOTLIB:
X                /* Y-cable for Tektornix? */
X      p1x = 0;    p1y = 0;
X      p2x = 5450;    p2y = 4095;
X                /* 5% room for labels */
X      t1 = p1x + 0.07*(p2x-p1x); /* y-label */
X      t2 = p2x - 0.02*(p2x-p1x);; /* righ margin */
X      xdevicelow  = t1 + xlow * (t2 - t1);
X      xdevicehigh = t1 + xhigh * (t2 - t1);
X      t1 = p1y + 0.10*(p2y-p1y); /* x-label */
X      t2 = p2y - 0.05*(p2y-p1y); /* name */
X      ydevicelow  = t1 + ylow * (t2 - t1);
X      ydevicehigh = t1 + yhigh * (t2 - t1);
X
X      openpl();
X      space(0,0,4096,4096);
X      erase();
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X
X    /* check coordinate position legality
X     */
X    if (xlow < -0.0001 || xhigh > 1.0001 ||
X    ylow < -0.0001 || yhigh > 1.0001) {
X    dr_errmsg("Illegal x/y low/high position.  Aborting.");
X    (void)dr_finish();
X    exit(1);
X    }
X    return;
X}
X
X
X/*
X * dr_xgrid  - draw x grid
X *
X * Draw grid lines parallel to y-axis.  When argument n == 0, makes
X * lines at every tic-mark, n == 1 makes one grid line between every
X * tic mark etc.  To leave grid lines away, just don't call this.
X */
Xvoid dr_xgrid(n)
Xint n;                /* number of grid lines */
X                /* between tic-marks */
X{
X    int i;
X    int x, y1, y2;        /* temps for grid line positions */
X    int nlines;            /* number of grid lines to draw */
X    char s[200];
X
X    /* check interval
X     */
X    if (n < 1) {
X    dr_errmsg("Invalid x-grid interval, using 1");
X    n = 1;
X    }
X
X    /* how many lines we have?
X     */
X    nlines = (int) (n*(xuserhigh-xuserlow)/xuserincr - 0.5);
X
X     switch (device)
X      {
X      case HP7475:
X      /*  select absolute moving
X       */
X      sends("PA;\n");
X#ifdef SLOW
X      if (slow) {
X      }
X#endif
X      /* draw desired number of lines, sawing up and down
X       */
X      y1 = yscale(yuserlow);
X      y2 = yscale(yuserhigh);
X      for (i=1; i<=nlines; i++) {
X                /* go up */
X          x = xscale(i*xuserincr/n+xuserlow);
X          (void)sprintf(s, "PU%d,%d;PD%d,%d;", x, y1, x, y2);
X          sends(s);
X                /* go down */
X          if (++i <= nlines) {
X          x = xscale(i*xuserincr/n+xuserlow);
X          (void)sprintf(s, "PU%d,%d;PD%d,%d;\n", x, y2, x, y1);
X          sends(s);
X          }
X      }
X
X      /* lift pen
X       */
X      sends("PU;\n");
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      y1 = yscale(yuserlow);
X      y2 = yscale(yuserhigh);
X      for (i=1; i<=nlines; i++) {
X          x = xscale(i*xuserincr/n+xuserlow);
X          n_line(x,y1, x,y2);
X      }
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      y1 = yscale(yuserlow);
X      y2 = yscale(yuserhigh);
X      for (i=1; i<=nlines; i++) {
X          x = xscale(i*xuserincr/n+xuserlow);
X          line(x,y1, x,y2);
X      }
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X}
X
X
X/*
X * dr_ygrid  - draw y grid
X *
X * See dr_xgrid().
X */
Xvoid dr_ygrid(n)
Xint n;
X{
X    int i;
X    int nlines;            /* number of lines */
X    char s[200];
X    int x1, x2, y;
X
X    /* check interval
X     */
X    if (n < 1) {
X    dr_errmsg("Invalid y-grid interval, using 1");
X    n = 1;
X    }
X
X    /* how many lines do we have?
X     */
X    nlines = (int) (n*(yuserhigh-yuserlow)/yuserincr - 0.5);
X
X    switch (device)
X      {
X      case HP7475:
X      /* select absolute moving
X       */
X      sends("PA;\n");
X
X      /* draw desired number of lines, sawing right and left
X       */
X      x1 = xscale(xuserlow);
X      x2 = xscale(xuserhigh);
X      for (i=1; i<=nlines; i++) {
X                /* go right */
X          y = yscale(yuserlow+i*yuserincr/n);
X          (void)sprintf(s, "PU%d,%d;PD%d,%d;", x1, y, x2, y);
X          sends(s);
X                /* go left */
X          if (++i <= nlines) {
X          y = yscale(yuserlow+i*yuserincr/n);
X          (void)sprintf(s, "PU%d,%d;PD%d,%d;", x2, y, x1, y);
X          sends(s);
X          }
X          sends("\n");
X      }
X
X      /* lift pen
X       */
X      sends("PU;\n");
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      x1 = xscale(xuserlow);
X      x2 = xscale(xuserhigh);
X      for (i=1; i<=nlines; i++) {
X          y = yscale(i*yuserincr/n+yuserlow);
X          n_line(x1,y, x2,y);
X      }
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      x1 = xscale(xuserlow);
X      x2 = xscale(xuserhigh);
X      for (i=1; i<=nlines; i++) {
X          y = yscale(i*yuserincr/n+yuserlow);
X          line(x1,y, x2,y);
X      }
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X}
X
X
X/*
X * dr_area_def  - establish x and y axis
X *
X * [xy]low    lower left corner
X * [xy]high   upper right corner
X * [xy]incr   increment of tic-marks
X * no_axes    flag, != 0 means don't draw axes
X *
X * Coordinate values will be printed at every tic-mark, so useful
X * values for incr is such that (high-low)/incr will be some small
X * integer.  With functions dr_xgrid() and dr_ygrid() it is possible to
X * draw grid lines also in between these tic-marcs.
X */
Xvoid dr_area_def(xlow, xincr, xhigh, ylow, yincr, yhigh, no_axes)
Xfloat xlow, xincr, xhigh;
Xfloat ylow, yincr, yhigh;
Xint no_axes;
X{
X    int i;
X    char units[30];        /* place for units */
X    char s[200];
X    float x, y;
X    int ix, iy, nx, ny;
X
X    /* set the globals
X     */
X    xuserlow = xlow;   xuserhigh = xhigh;   xuserincr = xincr;
X    yuserlow = ylow;   yuserhigh = yhigh;   yuserincr = yincr;
X    xfactor = (xdevicehigh - xdevicelow) / (xuserhigh - xuserlow);
X    yfactor = (ydevicehigh - ydevicelow) / (yuserhigh - yuserlow) ;
X
X    /* check legality,
X     * abort on error
X     */
X    if ((xincr == 0.0) ||
X    (xuserlow == xuserhigh) ||
X    ((xuserhigh-xuserlow)*xincr < 0.0)) {
X    dr_errmsg("x-axis defined improperly. Aborting.");
X    (void)dr_finish();
X    exit(1);
X    }
X    if ((yincr == 0.0) ||
X    (yuserlow == yuserhigh) ||
X    ((yuserhigh-yuserlow)*yincr < 0.0)) {
X    dr_errmsg("y-axis defined improperly. Aborting.");
X    (void)dr_finish();
X    exit(1);
X    }
X    if (no_axes)
X      return;
X
X    switch (device)
X      {
X      case HP7475:
X
X      /* draw box
X       * --------
X       */
X
X      /* select absolute plotting mode,
X       * lift pen,
X       * go to lower left corner
X       */
X      (void)sprintf(s,"PA;PU%d,%d;\n",
X            xscale(xlow), yscale(ylow));
X      sends(s);
X
X      /* draw x-axis with x-tick marks
X       */
X
X      (void)sprintf(s, "TL%.1f,%.1f;\n", 0.5*xrelarea, 0.5*yrelarea);
X      sends(s);
X      for (nx=0, x = xlow;
X           x < xhigh + xincr/100.0;
X           x += xincr, nx++) {
X          (void)sprintf(s, "PD%d,%d;XT;\n", xscale(x), yscale(ylow));
X          sends(s);
X      }
X      (void)sprintf(s, "PD%d,%d;\n", xscale(xhigh), yscale(ylow));
X      sends(s);
X
X      /* draw y-axis with y-tick marks
X       */
X      (void)sprintf(s, "PU%d,%d;\n",
X            xscale(xlow), yscale(ylow));
X      sends(s);
X      for (ny=0, y = ylow;
X           y < yhigh + yincr/100.0;
X           y += yincr, ny++) {
X          (void)sprintf(s, "PD%d,%d;YT;\n", xscale(xlow), yscale(y));
X          sends(s);
X      }
X      (void)sprintf(s, "PD%d,%d;\n", xscale(xlow), yscale(yhigh));
X      sends(s);
X
X      /* outer borders of the box
X       */
X      (void)sprintf(s, "PU%d,%d;PD%d,%d,%d,%d;PU\n",
X            xscale(xlow), yscale(yhigh),
X            xscale(xhigh), yscale(yhigh),
X            xscale(xhigh), yscale(ylow));
X      sends(s);
X      break;
X#ifdef GL_LIB
X#define GL_TIC_SIZE 200
X      case SCREEN:
X      iy = yscale(ylow);
X      n_movepen(xscale(xlow), iy);
X
X      for (nx=0, x = xlow;
X           x < xhigh + xincr/100.0;
X           x += xincr, nx++) {
X          ix = xscale(x);
X          n_draw(ix, iy);
X          n_line(ix, iy+GL_TIC_SIZE, ix, iy-GL_TIC_SIZE);
X          n_movepen(ix, iy);
X      }
X      n_draw(xscale(xhigh), iy);
X
X      /* draw y-axis with y-tick marks
X       */
X      ix = xscale(xlow);
X      n_movepen(ix, yscale(ylow));
X      for (ny=0, y = ylow;
X           y < yhigh + yincr/100.0;
X           y += yincr, ny++) {
X          iy = yscale(y);
X          n_draw(ix, iy);
X          n_line(ix+GL_TIC_SIZE, iy, ix-GL_TIC_SIZE, iy);
X          n_movepen(ix, iy);
X      }
X      n_draw(ix, yscale(yhigh));
X
X      /* outer borders of the box
X       */
X      n_movepen(xscale(xlow), yscale(yhigh));
X      n_draw(xscale(xhigh), yscale(yhigh));
X      n_draw(xscale(xhigh), yscale(ylow));
X      break;
X#endif
X#ifdef PLOT
X#define PLOT_TIC_SIZE 40
X      case PLOTLIB:
X      iy = yscale(ylow);
X      move(xscale(xlow), iy);
X
X      for (nx=0, x = xlow;
X           x < xhigh + xincr/100.0;
X           x += xincr, nx++) {
X          ix = xscale(x);
X          cont(ix, iy);
X          line(ix, iy+PLOT_TIC_SIZE, ix, iy-PLOT_TIC_SIZE);
X          move(ix, iy);
X      }
X      cont(xscale(xhigh), iy);
X
X      /* draw y-axis with y-tick marks
X       */
X      ix = xscale(xlow);
X      move(ix, yscale(ylow));
X      for (ny=0, y = ylow;
X           y < yhigh + yincr/100.0;
X           y += yincr, ny++) {
X          iy = yscale(y);
X          cont(ix, iy);
X          line(ix+PLOT_TIC_SIZE, iy, ix-PLOT_TIC_SIZE, iy);
X          move(ix, iy);
X      }
X      cont(ix, yscale(yhigh));
X
X      /* outer borders of the box
X       */
X      move(xscale(xlow), yscale(yhigh));
X      cont(xscale(xhigh), yscale(yhigh));
X      cont(xscale(xhigh), yscale(ylow));
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X
X
X    /* put y-units
X     */
X    for (i=ny-1; i>=0; i--) {
X     (void)sprintf(units, "%g", ylow + i * yincr);
X     dr_put_text(xlow, ylow+i*yincr,
X            -strlen(units)/2.0, -1.3, -M_PI_2, 1.0, units);
X    }
X
X    /* put x-units
X     */
X    for (i=0; i<nx; i++) {
X    (void)sprintf(units, "%g", xlow + i * xincr);
X    dr_put_text(xlow+i*xincr, ylow,
X            -strlen(units)/2.0, -1.3, 0.0, 1.0, units);
X    }
X    return;
X}
X
X/*
X * linetypes
X *
X */
X#ifdef GL_LIB
Xstatic long gl_line_types[] = {    /* 76543210765432107654321076543210 */
X    0xf0000000L,        /* ====............................ */
X    0xffff0000L,        /* ================................ */
X    0xffffff00L,        /* ========================........ */
X    0xffffff18L,        /* ========================...==... */
X    0xfffff0f0L,        /* ====================....====.... */
X    0xfffff39cL,        /* ====================..===..===.. */
X    0xffffffffL,        /* ================================ */
X    0x66666666L,        /* ==..==..==..==..==..==..==..==.. */
X    0xf8f8f8f8L,        /* =====...=====...=====...=====... */
X    0xffffff00L,        /* ========================........ */
X    0xfc30fc30L };        /* ======....==....======....==.... */
X#define DEF_GL_TYPE 7        /* this will produce continous line */
X#endif /* GL_PLOT */
X#ifdef PLOT
Xstatic char *plot_line_types[] = {
X    "dotted",            /* 1 */
X    "longdashed",        /* 2 */
X    "shortdashed",        /* 3 */
X    "dotdashed",        /* 4 */
X    "solid",            /* 5 */
X    "solid",            /* 6 */
X    "solid" };            /* 7 */
X#define DEF_PLOT_TYPE 7        /* this will produce continous line */
X#endif /* PLOT */
X
Xint dr_set_ltype(linetype)
Xint linetype;
X{
X    char s[40];
X    int old_val = old_linetype;
X
X    switch (device)
X      {
X      case HP7475:
X      if (linetype < 0 || linetype > 6) {
X          linetype = 7;
X          sends("LT;");
X      }
X      else {
X          (void)sprintf(s, "LT%d;", linetype);
X          sends(s);
X      }
X      break;
X#ifdef GL_LIB
X      case SCREEN:
X      if (linetype < 1 ||
X          linetype > (sizeof(gl_line_types) / sizeof(gl_line_types[0]))) {
X           sprintf (s, "Invalid linetype: %d", linetype);
X           dr_errmsg(s);
X           linetype = DEF_GL_TYPE;
X       }
X      g_style(gl_line_types[linetype-1]);
X      break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X      if (linetype < 1 ||
X          linetype > (sizeof(plot_line_types) /
X              sizeof(plot_line_types[0]))) {
X           sprintf (s, "Invalid linetype: %d", linetype);
X           dr_errmsg(s);
X           linetype = DEF_PLOT_TYPE;
X       }
X      linemod(plot_line_types[linetype-1]);
X      break;
X#endif
X      default:
X      invdev(); break;
X      }
X    old_linetype = linetype;
X    return old_val;
X}
X
X/*
X * xscale  - scaling in horizontal directon
X *
X * convert user x-coordinate to physical device x-coordinate
X */
Xstatic int xscale(x)
Xfloat x;
X{
X    int temp;
X    char s[60];
X
X    /* calculate it
X     */
X    temp = (int) ((x - xuserlow) * xfactor + xdevicelow + 0.5);
X    /* check, if it is good
X     */
X    if ((temp >= -32768) && (temp <= 32767))
X      return (temp);
X
X    /* when not good make report,
X     * return something
X     */
X    (void)sprintf(s, "x-value out of range: %f", x);
X    dr_warnmsg(s);
X    return (int) (xdevicelow + 0.5);
X}
X
X
X/*
X * yscale  - scaling in vertical directon
X *
X * convert user y-coordinate to physical device y-coordinate
X */
Xstatic int yscale(y)
Xfloat y;
X{
X    int temp;
X    char s[60];
X
X    /* calculate it
X     */
X    temp = (int) ((y - yuserlow) * yfactor + ydevicelow + 0.5);
X
X    /* check, if it is good
X     */
X    if ((temp >= -32768) && (temp <= 32767))
X      return (temp);
X
X    /* when not good make report,
X     * return something
X     */
X    (void)sprintf(s, "y-value out of range: %f", y);
X    dr_warnmsg(s);
X    return (int) (ydevicelow + 0.5);
X}
X
Xstatic void invdev()
X{
X    dr_errmsg("Invalid plotting device.");
X    (void)dr_finish();
X    exit(1);
X}
X
X
X/*
X * dr_errmsg  - send a string to error output
X *
X */
Xvoid dr_errmsg(s)
Xchar *s;
X{
X    char *p = s + strlen(s);
X    if (p != s) p--;
X    while ((p != s) && (*p == '\n')) /* strip the newlines */
X      *(p--) = '\0';
X    fprintf(errfilefp, "dr_lib error: %s\n", s);
X    errcount++;
X    return;
X}
X
X
X/*
X * warnmsg  - send a string to error output
X *
X */
Xvoid dr_warnmsg(s)
Xchar *s;
X{
X    char *p = s + strlen(s);
X    if (p != s) p--;
X    while ((p != s) && (*p == '\n')) *(p--) = '\0';
X    fprintf(errfilefp, "Warning: %s\n", s);
X    warncount++;
X    return;
X}
X
X
X/*
X * sends  - send something to plotting device
X *
X */
Xstatic void sends(s)
Xchar *s;
X{
X    (void)printf("%s", s);
X}
X
X#ifdef ASK
X/*
X * takes  - take a string from plotter (to read plotter coordinates)
X *
X */
X#include <signal.h>
X#include <setjmp.h>
X
Xjmp_buf sjbuf;
Xint (*savAlrm)();        /* save incomming alarm function */
X
Xstatic int timerh()        /* timer interrupt handler */
X{
X    longjmp(sjbuf, 1);
X}
X
Xstatic void takes(s)
Xchar *s;
X{
X#ifdef YCABLE
X    savAlrm = signal(SIGALRM,timerh);
X    alarm(2);            /* give time to wake up */
X    if (setjmp(sjbuf)) {
X    alarm(0);
X    signal(SIGALRM,savAlrm); /* cancel timer */
X    (void)dr_finish();
X    (void)fprintf(stderr, "\n\7Timeout reading stdin.  Aborting.\n");
X    exit(1);
X    }
X    gets(s);
X    alarm(0);
X    signal(SIGALRM,savAlrm);    /* cancel timer */
X#else
X    /* do it yorself */
X#endif
X}
X#endif /* ASK */
END_OF_FILE
if test 31312 -ne `wc -c <'dr_lib.c'`; then
    echo shar: \"'dr_lib.c'\" unpacked with wrong size!
fi
# end of 'dr_lib.c'
fi
if test -f 'draw.c' -a "${1}" != "-c" ; then
  echo shar: Will not clobber existing file \"'draw.c'\"
else
echo shar: Extracting \"'draw.c'\" \(17080 characters\)
sed "s/^X//" >'draw.c' <<'END_OF_FILE'
X/*
X
XCopyright (C) 1988, 1989 by Juha Takala, [e-mail]
X
X     This program is free software; you can redistribute it and/or modify
X     it under the terms of the GNU General Public License as published by
X     the Free Software Foundation; version 1.
X
X     This program is distributed in the hope that it will be useful,
X     but WITHOUT ANY WARRANTY; without even the implied warranty of
X     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
X     the file "License for more details
X
X     You should have received a copy of the GNU General Public License
X     along with this program; if not, write to the Free Software
X     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
X
X
X */
X
X/* Description:
X *  Non-interactive user interface program to provide interface between
X *  device manipulation subroutines and user supplied command files.
X *  This program will accept all kinds of command line arguments and
X *  pseudo comments mixed in between the data lines.  See the man page.
X *
X *  The good thing about this program is that it will accept infinite
X *  amount of input data points; they are not stored, but processed one
X *  at a time.
X */
Xstatic char *RCS_id =
X  "$Id: draw.c,v 2.7 89/12/08 14:00:39 jta Exp Locker: jta $";
X
X#include <stdio.h>
X#ifdef MSDOS
X#include <stdlib.h>
X#endif /* MSDOS */
X#include <string.h>
X#include <sys/types.h>
X#include <time.h>
X#include <math.h>
X#include <signal.h>
X#ifndef MSDOS
X#include <sys/inode.h>
X#include <sys/file.h>
X#endif
X#include <string.h>
X#include "dr_lib.h"
X
Xdouble strtod();
X
Xchar *myname;
Xfloat xlowl = 0.0;        /* relative plotting area on paper */
Xfloat ylowl = 0.0;
Xfloat xhighr= 1.0;
Xfloat yhighr = 1.0;
X
Xfloat xlow = -1.0;        /* x axis range */
Xfloat xincr = 0.5;
Xfloat xhigh = 1.0;
Xint xgrid_n = 0;
Xchar *xlabel_text = NULL;
X
Xfloat ylow = -1.0;        /* y axis range */
Xfloat yincr = 0.5;
Xfloat yhigh = 1.0;
Xint ygrid_n = 0;
Xchar *ylabel_text = NULL;
X
Xstatic float x, y;        /* current point, needs to be */
X                /* remembered, eg. when relative text is */
X                /* inserted by #$texthere */
Xchar *name_text = NULL;        /* name of picture */
Xint pennum = 1;
Xint penup = 1;
Xint no_axes = 0;
Xint init_done = 0;
Xint dateflag = 0;
Xint linetype = 7;        /* default for solid linetype */
Xfloat text_angle = 0.0;        /* radians */
Xfloat text_size = 1.0;        /* text size factor */
Xfloat text_shape = 1.0;        /* text shape factor */
Xint cset = 0;            /* char set, default is US Ascii */
Xint continuous = 1;        /* continuous line / discrete points */
X                /* 0 = discontinous */
Xchar marker = '*';        /* default marker char in discontinous mode */
Xchar informat[200] = "%f %f";    /* default input data line format */
Xint reversed = 0;        /* != 0 means input is y before x  */
Xint absolute = 0;        /* != 0 when 'absolute' coordinates are */
X                /* to be processed */
Xfloat xabs_origin = 0.0;    /* origin of 'absolute' mode data points */
Xfloat yabs_origin = 0.0;
X
Xint device = HP7475 | SIZE_A4;
X
Xstatic int setmode();
Xstatic int labelpic();
Xstatic void doit();
Xstatic void setdevice();
Xextern void exit();
X
Xstatic void (*old_int_handler)();
Xstatic void new_int_handler();
Xstatic void fatal();
Xstatic void usage();
X
X#define PUT_TEXT(x,y,xofs,yofs,text) \
X      dr_put_text((x), (y), (xofs), (yofs), text_angle, 1.0, (text))
X#define PUT_X_TEXT(x,y,text) \
X      dr_put_text((x), (y), -(float)strlen(text), -2.0, 0.0, 1.4, (text))
X#define PUT_Y_TEXT(x,y,text) \
X      dr_put_text((x), (y), 0.0, -2.0, -M_PI_2, 1.4, (text))
X#define PUT_N_TEXT(x,y,text) \
X      dr_put_text((x), (y), -(float)strlen(text)/2.0, 0.4, 0.0, 2.0, (text))
X
Xmain(argc, argv)
Xint argc;
Xchar *argv[];
X{
X    int c;
X    extern char *optarg;
X    extern int optind, opterr;
X    char *opts = "a:bc:df:g:hl:m:n:p:s:S:t:x:X:y:Y:v";
X
X    FILE *fp;
X
X    myname = *argv;
X    setdevice(getenv("DRAWDEV"));
X
X
X    while ((c = getopt(argc, argv, opts)) != -1) {
X    switch (c) {
X    case 'a':
X        if (sscanf(optarg, "%f,%f,%f,%f",
X               &xlowl, &ylowl, &xhighr, &yhighr) != 4)
X          fatal("invalid plotting area definition: ", optarg);
X        break;
X    case 'b': no_axes++; break;
X    case 'c': cset = atoi(optarg); break;
X    case 'd': dateflag++; break;
X    case 'f':
X        /* input lines are : "-- -- y -- -- x -- --"
X         * (dashes indicate fields, that we are not interested in)
X         * and we want to plot y=f(x), so we specify option
X         * "-f 6,3" (fields are numbered from 1...)
X         */
X        {
X        int xcol, ycol;
X        if (sscanf(optarg, "%d,%d", &xcol, &ycol) != 2)
X          fatal("invalid input format definition: %s", optarg);
X        if (xcol == ycol)
X          fatal("invalid input format definition: %s", optarg);
X        if (xcol > ycol) {
X            int temp = xcol;
X            reversed++;
X            xcol = ycol;
X            ycol = temp;
X        }
X        ycol -= xcol;
X        informat[0] = '\0';
X        while (--xcol)
X          (void)strcat(informat, "%*f ");
X        (void)strcat(informat, "%f ");
X        while (--ycol)
X          (void)strcat(informat, "%*f ");
X        (void)strcat(informat, "%f");
X        }
X        break;
X    case 'g':
X        if (sscanf(optarg, "%d,%d", &xgrid_n, &ygrid_n) != 2)
X          fatal("invalid grid spec.: ", optarg);
X        break;
X    case 'h': case '?':
X        usage(); break;
X    case 'l':
X        if (sscanf(optarg, "%d", &linetype) != 1)
X          fatal("invalid line type:", optarg);
X        break;        /* draw axes before setting it! */
X    case 'm':
X        continuous = 0;
X        if (sscanf(optarg, "%1s", &marker) != 1)
X          fatal ("no marker found.", "\0");
X        break;
X    case 'n':
X        name_text = optarg; break;
X    case 'p':
X        if (sscanf(optarg, "%d", &pennum) != 1)
X          fatal("invalid pen:", optarg);
X        if ((pennum > 6) || (pennum < 1))
X          pennum = 1;
X        break;
X    case 's':
X        if (sscanf(optarg, "%f", &text_size) != 1)
X          fatal("invalid text size: ", optarg);
X        break;
X    case 'S':
X        if (sscanf(optarg, "%f", &text_shape) != 1)
X          fatal("invalid text shape: ", optarg);
X        break;
X    case 't': setdevice(optarg); break;
X    case 'x':
X        if (sscanf(optarg, "%f,%f,%f",
X               &xlow, &xincr, &xhigh) != 3)
X          fatal("invalid x range: ", optarg);
X        break;
X    case 'X':
X        xlabel_text = optarg; break;
X    case 'y':
X        if (sscanf(optarg, "%f,%f,%f",
X               &ylow, &yincr, &yhigh) != 3)
X          fatal("invalid y range: ", optarg);
X        break;
X    case 'Y':
X        ylabel_text = optarg; break;
X    case 'v':
X        fprintf(stderr, "%s\n", RCS_id);
X        exit (0);
X    default:
X        exit(1);
X    }
X    }
X    argc -= optind;
X    argv += optind;
X
X    old_int_handler = signal (SIGINT, new_int_handler);
X
X    /*
X     * Options have been taken, start the work.  Initializations will be
X     * done just before first data point is used...  This is jippo to
X     * make it possible to give coordinate specifications from inside
X     * the data....
X     */
X    if (argc == 0)        /* no files, take points from stdin */
X      doit(stdin);
X    else {            /* take points from all files */
X    for ( ; argc; argc--, argv++) {
X        if (NULL == (fp = fopen(*argv, "r"))) {
X        perror(*argv);
X        continue;
X        }
X        doit(fp);
X        fclose(fp);
X    }
X    }
X    if (init_done)
X      return(dr_finish(1) ? 1 : 0); /* 1 == error(s) */
X    else
X      return 2;            /* nothing was done! */
X}
X
X
X/* gl_plot will set int handler at its initialization time, and restore
X   what it saw there previously, and finally it will
X   kill(SIGINT,getpid()).  This is what we are going to do, too, so that
X   OUR caller can clean ITS data stuctures etc.
X*/
Xstatic void new_int_handler() {
X    (void)dr_finish(0);
X    signal (SIGINT, old_int_handler);
X    kill (getpid(), SIGINT);    /* zap my self */
X    sleep(1);
X    exit(3);            /* if that is not enough */
X}
X
X
Xstatic void do_initializations()
X{
X    init_done++;        /* remember it is done */
X    dr_start(xlowl, ylowl, xhighr, yhighr, device);
X    if (cset != 0) dr_set_cset(cset);
X    dr_set_size(text_size);
X    dr_set_shape(text_shape);
X    dr_set_pen(pennum);
X    dr_area_def(xlow, xincr, xhigh, ylow, yincr, yhigh, no_axes);
X    if (!no_axes) {
X    if (xgrid_n) dr_xgrid(xgrid_n);
X    if (ygrid_n) dr_ygrid(ygrid_n);
X    }
X    if (xlabel_text) PUT_X_TEXT(xhigh, ylow, xlabel_text);
X    if (ylabel_text) PUT_Y_TEXT(xlow, yhigh, ylabel_text);
X    if (name_text) PUT_N_TEXT((xhigh+xlow)/2.0, yhigh, name_text);
X    if (dateflag) putdate(dateflag);
X    dr_set_ltype(linetype);
X}
X
Xstatic void setdevice(s)
Xchar *s;
X{
X    if (s == NULL) return;
X    if      (!strcmp(s,"a4")) device = (HP7475 | SIZE_A4);
X    else if (!strcmp(s,"a3")) device = (HP7475 | SIZE_A3);
X    else if (!strcmp(s,"a"))  device = (HP7475 | SIZE_A);
X    else if (!strcmp(s,"b"))  device = (HP7475 | SIZE_B);
X#ifdef GL_LIB
X    else if (!strcmp(s,"ega")) device = SCREEN;
X#endif
X#ifdef PLOT
X    else if (!strcmp(s,"plot")) device = PLOTLIB;
X#endif
X#ifdef TEK
X    else if (!strcmp(s,"tek"))
X      fatal("Sorry,", "Tek4010 not yet.");
X#endif
X    else fatal("Unsupported plotting device", s);
X}
X
Xstatic void doit(fp)
XFILE *fp;
X{
X#define LINESIZE 512
X    char buf[LINESIZE+1], *p;
X    char emsg[200];
X    int lineno = 0;
X
X    while (1) {
X    if (NULL == fgets(buf, LINESIZE, fp)) /* EOF */
X      return;
X
X    lineno++;
X    for (p=buf; *p; p++)    /* skip white space */
X      if (strchr(" \t",*p) == NULL)
X        break;
X    if (*p == '\n')
X      continue;
X
X    if (*p == '#') {
X        if (setmode(p+1)) {
X        sprintf(emsg, "Invalid pseudo comment line (%d): %s",
X            lineno, buf);
X        if (init_done)
X          dr_errmsg(emsg);
X        else
X          write(2, emsg, strlen(emsg));
X        }
X        continue;
X    }
X                /* now we expect data lines, */
X                /* initializations must be done by now */
X    if (!init_done) do_initializations();
X
X    if (sscanf(p, informat, &x, &y) != 2) {
X        sprintf(emsg, "Invalid data line (%d): %s", lineno, buf);
X        dr_errmsg(emsg);
X        continue;
X    }
X    if (reversed) {
X        float temp = x;
X        x = y;
X        y = temp;
X    }
X    if (absolute) {
X        x = xlow + (xabs_origin + x) * (xhigh - xlow);
X        y = ylow + (yabs_origin + y) * (yhigh - ylow);
X    }
X    /* this needs to be fixed: reversed directions not handled! */
X    if (x < xlow || x > xhigh || y < ylow || y > yhigh) {
X        sprintf (emsg, "point out of range at line (%d): %s",
X             lineno, buf);
X        dr_warnmsg(emsg);
X    }
X    if (continuous) {
X        dr_goto(x, y, penup);
X        penup = 0;
X    }
X    else {
X        dr_put_mark(x, y, marker);
X    }
X    }
X}
X
Xstatic void remember(to, from)    /* copy string to safe place, put its */
X                /* addr into `*to' */
Xchar **to, *from;
X{
X#ifndef MSDOS
X    extern char *malloc();
X#endif
X    *to = malloc(strlen(from)+1);
X    strcpy(*to, from);
X}
X
X                /* handle '#<something>' lines */
Xstatic int setmode(s)        /* return 0 on success, != when trouble */
Xchar *s;
X{
X    for (; *s; s++)        /* skip white space */
X      if (strchr(" \t",*s) == NULL)
X    break;
X    if (*s++ != '$')        /* only '$' causes action */
X      return 0;
X
X    for (; *s; s++)        /* skip white space */
X      if (strchr(" \t",*s) == NULL)
X    break;
X
X    if (!strncmp (s, "absolute", 8)) {
X    absolute = 1;
X    if (sscanf(s+8, "%f,%f", &xabs_origin, &yabs_origin) == 2)
X      return 0;
X    xabs_origin = 0.0;
X    yabs_origin = 0.0;
X    return 1;
X    }
X
X    if (!strncmp (s, "endabsolute", 11)) {
X    absolute = 0;
X    return 0;
X    }
X
X    if (!strncmp (s,"penup",5)) {
X    penup = 1;
X    return 0;
X    }
X
X
X    if (!strncmp(s, "include", 7)) {
X    FILE *fp;
X    char *p;
X
X    for (s+=8; *s; s++)    /* skip whitespace */
X      if (strchr(" \t",*s) == NULL)
X         break;
X    for (p=s; *p; p++)
X       if (strchr(" \t\n",*s) != NULL)
X         break;
X     *(p-1) = '\0';        /* put terminator */
X
X     if (NULL == (fp = fopen(s, "r"))) {
X         perror(p);
X         return 1;
X     }
X     doit(fp);        /* recursive call */
X     fclose(fp);
X     return 0;
X    }
X
X    if (!strncmp (s, "cset", 4)) {
X    int csetnum = 0;
X    if (sscanf(s+4, "%d", &csetnum) != 1)
X      return 1;
X    if (!init_done) do_initializations();
X    dr_set_cset(csetnum);
X    return 0;
X    }
X
X    if (!strncmp (s, "pen", 3)) {
X    if (sscanf(s+3, "%d", &pennum) != 1)
X      return 1;
X    if (!init_done) do_initializations();
X    dr_set_pen(pennum);
X    return 0;
X    }
X
X    if (!strncmp (s, "ltype", 5)) {
X    int linetype = 1;
X    if (sscanf(s+5, "%d", &linetype) != 1)
X      return 1;
X    if (!init_done) do_initializations();
X    dr_set_ltype(linetype);
X    return 0;
X    }
X
X    if (!strncmp(s, "mark", 4)) {
X    s += 4;
X    while (1) {
X        marker = *s++;
X        if (strchr(" \t", marker) == NULL)
X          break;
X    }
X    continuous = (strchr (" \n\r\t", marker) != NULL);
X    return 0;
X    }
X
X    if (!strncmp(s, "date",   4)) {
X    if (dateflag)
X      return 1;
X    dateflag = 1;
X    if (init_done)
X      putdate(1);
X    return 0;
X    }
X    if (!strncmp(s, "time",   4)) {
X    if (dateflag)
X      return 1;
X    dateflag = 2;
X    if (init_done)
X      putdate(2);
X    return 0;
X    }
X
X    if (!strncmp(s, "xlabel", 6)) {
X    if (xlabel_text) return 1;
X    s += 7;            /* skip one space/separator */
X    s[strlen(s)-1] = '\0';    /* drop newline */
X    if (!init_done) {
X        remember(&xlabel_text, s);
X        return 0;
X    }
X    PUT_X_TEXT(xhigh, ylow, s);
X    xlabel_text = myname;    /* some char-pointer that has non NULL */
X                /* value, so that we remember this */
X    penup = 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "ylabel", 6)) {
X    if (ylabel_text) return 1;
X    s += 7;
X    s[strlen(s)-1] = '\0';
X    if (!init_done) {
X        remember(&ylabel_text, s);
X        return 0;
X    }
X    PUT_Y_TEXT(xlow, yhigh, s);
X    ylabel_text = myname;
X    penup = 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "name", 4)) {
X    if (name_text) return 1;
X    s += 5;            /* skip the space */
X    s[strlen(s)-1] = '\0';    /* drop newline */
X    if (!init_done) {
X        remember(&name_text, s);
X        return 0;
X    }
X    PUT_N_TEXT((xhigh+xlow)/2.0, yhigh, s);
X    name_text = myname;
X    penup = 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "xrange", 6)) {
X    if (sscanf(s+6, "%f,%f,%f", &xlow, &xincr, &xhigh) != 3)
X      return 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "yrange", 6)) {
X    if (sscanf(s+6, "%f,%f,%f", &ylow, &yincr, &yhigh) != 3)
X      return 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "area", 4)) {
X    if (sscanf(s+4, "%f,%f,%f,%f",
X           &xlowl, &ylowl, &xhighr, &yhighr) != 4)
X      return 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "grid", 4)) {
X    if (sscanf(s+4, "%d,%d", &xgrid_n, &ygrid_n) != 2)
X      return 1;
X    return 0;
X    }
X
X    if (!strncmp(s, "noaxes", 6)) {
X    no_axes++;
X    return 0;
X    }
X
X    if (!strncmp(s, "init", 4)) {
X    if (!init_done) do_initializations();
X    return 0;
X    }
X
X    if (!strncmp(s, "size", 4)) {
X    if (sscanf(s+4, "%f", &text_size) != 1)
X      return 1;
X    if (init_done)
X      dr_set_size(text_size);
X    return 0;
X    }
X
X    if (!strncmp(s, "shape", 5)) {
X    if (sscanf(s+5, "%f", &text_shape) != 1)
X      return 1;
X    if (init_done)
X      dr_set_shape(text_shape);
X    return 0;
X    }
X
X    if (!strncmp(s, "direction", 9)) {
X    float xx;
X
X    if (sscanf(s+9, "%f", &xx) != 1)
X      return 1;
X    text_angle = xx;
X    return 0;
X    }
X
X    if (!strncmp(s, "text",   4)) return (labelpic(s+4, 0, 0));
X    if (!strncmp(s, "alabel", 6)) return (labelpic(s+6, 0, 1));
X    if (!strncmp(s, "rlabel", 6)) return (labelpic(s+6, 1, 1));
X
X                /* other options are errors */
X    return 1;
X}
X
Xstatic int labelpic(s, rel, conn)
Xchar *s;
Xint rel, conn;
X{
X    float xx, yy;
X    char *s1;
X    int old_linetype;
X
X    if (!init_done) do_initializations();
X
X    old_linetype = dr_set_ltype(7); /* set solid line */
X    xx = strtod(s, &s1);
X    if (s1 == s) return 1;
X    s = s1+1;
X    yy = strtod(s, &s1);
X    if (s1 == s) return 1;
X    s1++;            /* skip the separator */
X    s1[strlen(s1)-1] = '\0';
X    if (rel) {
X    xx += x;
X    yy += y;
X    }
X    PUT_TEXT(xx, yy, 0.0, 0.15, s1);
X    if (conn) {
X    dr_goto(xx, yy, 0);
X    dr_goto(x, y, 0);
X    dr_draw_circle();
X    }
X    else
X      penup = 1;
X    (void) dr_set_ltype(old_linetype);
X    return 0;
X}
X
Xint putdate(flag)
Xint flag;
X{
X    struct tm *timp;
X    long clk;
X    char buf1[40], buf2[20];
X
X    if (time(&clk) < 0) {
X    (void) sprintf (buf1, "time() failed.");
X    (void) dr_errmsg(buf1);
X    return 1;
X    }
X    timp = localtime(&clk);
X    sprintf(buf1, "%d/%d/%d",
X        timp->tm_mday, timp->tm_mon+1, timp->tm_year + 1900);
X    if (flag > 1)
X      sprintf(buf2, " %02d:%02d:%02d",
X          timp->tm_hour, timp->tm_min, timp->tm_sec);
X    else
X      buf2[0] = '\0';
X    strcat(buf1, buf2);
X    PUT_TEXT(xhigh, yhigh, -(float)strlen(buf1), 0.7, buf1);
X    penup = 1;
X    return 0;
X}
X
X
Xvoid fatal(s1, s2)
Xchar *s1, *s2;
X{
X    fprintf (stderr,"%s: Fatal error: %s %s\n", myname, s1, s2);
X    exit(1);
X}
X
Xstatic char *utbl[] = {
X    "-a xlo,ylo,xhi,yhi\tspecify relative plot area on paper",
X    "-b\t\t\tdon't draw axes",
X    "-c num\t\t\tchoose character set `num'",
X    "-d\t\t\tinclude date (time) stamp",
X    "-f c1,c2\t\tuse these columns of input, higher is for y",
X    "-g xgrid,ygrid\t\tspecify x&y grids",
X    "-h\t\t\tgive this stuff",
X    "-l\t\t\tlinetype",
X    "-m marker\t\tspecify marker & discontinous mode",
X    "-n text\t\tname the picture",
X    "-p num\t\t\tspecify pen number",
X#ifdef GL_LIB
X    "-t target\t\tspecify device: `ega', `a3', `a4' etc. Defaults from env.",
X#else
X    "-t target\t\tspecify device: `a3', `a4' etc. Defaults from env.",
X#endif
X    "-x|y low,incr,high\tspecify axis range",
X    "-X|Y text\t\taxis label",
X    "-v\t\t\tprint version ID",
X    "Most (all?) of these may also be given inside the data, look for manual",
X    NULL
X};
X
Xvoid usage()
X{
X    register char **p;
X    printf ("usage: \"%s [opts] [file ...]\"\n", myname);
X    for (p=utbl; *p; p++)
X      (void)printf (" %s\n", *p);
X    exit(0);
X}
END_OF_FILE
if test 17080 -ne `wc -c <'draw.c'`; then
    echo shar: \"'draw.c'\" unpacked with wrong size!
fi
# end of 'draw.c'
fi
echo shar: End of archive 2 \(of 2\).
cp /dev/null ark2isdone
MISSING=""
for I in 1 2 ; do
    if test ! -f ark${I}isdone ; then
    MISSING="${MISSING} ${I}"
    fi
done
if test "${MISSING}" = "" ; then
    echo You have unpacked both archives.
    rm -f ark[1-9]isdone
else
    echo You still need to unpack the following archives:
    echo "        " ${MISSING}
fi
##  End of shell archive.
exit 0